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SHORT ARC OPTICAL SURVEY TECHNIQUES 

by 

J. Berbert and F. Loveless 
ABSTRACT 

The effect of the gravity parameter, p, the choice of fixed origin station, 
the local survey of the fixed origin station, and the choice of initial datum on 
the results of short arc satellite survey adjustments are investigated using 
GEOS-1 MOTS optical observations from 13 stations. It is concluded that each 
of these parameters has an effect on derived network scale on the order of 
2 x 10 6 for the nominal variations used in this study. Consequently, best 
available values should be used. A particular solution using what we thought to 
be the best available values for these parameters is recommended. 
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SHORT ARC OPTICAL SURVEY TECHNIQUES 


1. INTRODUCTION 

The purpose of this report is to present the results of studies to determine 
the effect of various parameters on the short arc recovery of survey for 13 
GEOS-1 MOTS optical stations using real data. The parameters investigated 
are 

(1) The error in the gravity parameter, fi, 

(2) The choice of fixed origin station, 

(3) The error in the assumed position of the fixed origin station, 

(4) The choice of datum to which the initial station positions are referred. 

This report documents and further details the results of geodetic studies for 
which preliminary results were given in References 1, 2, 3, and 4. 

Survey solutions are obtained for different values of each of the investigated 
parameters to illustrate the sensitivity of the solution to changes in these pa- 
rameters. A best set of parameters is chosen and the resulting survey solution 
is compared with the earlier results given by Brown (Reference 1), where a dif- 
ferent set of parameters and a somewhat larger network were used. 

The goal of this work is to develop the techniques for short arc survey ad- 
justment using MOTS camera GEOS satellite tracking data, and if possible, to 
reduce the uncertainty in the conventional surveys for these stations. The un- 
certainty in the relative positions of the MOTS stations represents a significant 
portion of the error budgets of the optical reference short arcs generated for 
use in the GEOS Observation Systems Intercomparison Investigation. 


2. THE OPTICAL TRACKING NETWORK 

The network includes 13 NASA MOTS camera stations as depicted in Figure 
1. In all our short arc solutions, the position of one MOTS camera, usually the 
centrally located camera at Columbia, Missouri, was held fixed. All other cam- 
era positions were allowed to freely adjust along with the six short arc orbital 
elements of each of the 38 GEOS-1 passes observed. The ground tracks of these 
selected passes and their relationship to the tracking stations are illustrated in 
Figure 2. As many as seven different well-spaced, 7-flash sequences were ob- 
served on some passes, but the typical number was four or five. In general, 
each camera was exercised to the maximum practical extent, so that as many 
as four separate 7 -flash sequences were sometimes obtained from a single cam- 
era on a single pass. 
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Figure 2. Ground Traces Relative to Tracking Network of Set of 38 Passes 
Carried in Short Arc Reduction 
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3. DATA REDUCTION 


The GEOS-1 MOTS camera observations were obtained on cards from the 
NSSDC at Goddard. The given UTC flash observation times were transformed 
to UT1 and 0.5 milliseconds were added to account for flash buildup time. The 
given right ascensions and declinations were corrected to allow for polar motion 
relative to the adopted IPMS mean pole of 1900-1905. 

Simultaneous multipass short arc survey adjustments were done with the 
Geodetic Data Adjustment Program (GDAP). This program was described by 
Lynn in Reference 6. Briefly, it is an orbit integration and generalized least 
squares parameter recovery program. The orbit integrations in GDAP are 
performed by means of power series expansions developed by Hartwell (Refer- 
ence 7). GDAP is capable of recovering an unlimited number of orbits simul- 
taneously with up to 64 error model or survey parameters. In these solutions, 
the recovered parameters were the East, North, and Up (AE, AN, AV) topocentric 
corrections to the initial estimates of the coordinates of the 12 adjusted stations, 
a total of 3 x 12 = 36 parameters. Two or three iterations of the solution of the 
normal equations are usually sufficient to minimize the sum of squares of the 
weighted obseiwation residuals and to achieve convergence to the final estimates 
of the survey adjustments. 

As mentioned by Brown (Reference 1), the errors in the adopted geopotential 
coefficients can be almost totally absorbed by the recovered orbital elements, 
so that short arc trajectory errors due to geopotential errors can be held to 
about one meter. Thus, the GDAP geopotential model, which includes only the 
zonal harmonics, J 2 through J 7 , is probably sufficient for these reductions. The 
zonal harmonic coefficients used are given below. They were taken from the 
NWL 5E-6 solutions (Reference 8) where they are associated with a gravity con- 
stant, n, of 3.986 0542 x 10 14 m 3 /sec 2 . 


Denormalized Zonal Harmonic Coefficients (xlO ' 6 ) 


J2 

J3 

J4 

J5 

J6 

J7 

1082.640 

-2.546 

-1.649 

-0.210 

0.640 

-0.333 


Survey solutions were obtained to determine the effect of variations in the 
parameters investigated, as indicated in Table 1 below. 


Tables 2, 3, and 4 list the initial survey coordinates in geodetic latitude, 
longitude, and geodetic height (i p, h), and in earth fixed geocentric cartesian 

coordinates (X, Y, Z) for the NAD-27, the SAO C-5, and the Mercury datums. 
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Table 1 


Parameter Variations for Each Solution 


Parameter 

Investigated 

Solution 

Gravity 
Parameter,/* 
(xio 14 m 3 /sec 2 ) 

Origin 

Datum 

Effect of error in /* 

1 

3.986032 

Columbia 

SAO-C5 


2 

3.986013 

Columbia 

SAO-C5 


3 

3.985994 

Columbia 

SAO-C5 


4 

3.985956 

Columbia 

SAO-C5 

Effect of Origin Location 

4 

3.985956 

Columbia 

SAO-C5 


5 

3.985956 

Rosman 

SAO-C5 


6 

3.985956 

Mojave 

SAO-C5 

Effect of error in Origin 

7 

3.986032 

Columbia(X o , Yo , Z o ) 

Mercury 


8 

3.986032 

Columbia (X o -50m) 

Mercury 


9 

3.986032 

Columbia (Y o -50m) 

Mercury 


10 

3.986032 

Columbia (Y o +50m) 

Mercury 


11 

3.986032 

Columbia (Z o -50m) 

Mercury 

Effect of initial Datum* 

2 

3.986013 

Columbia 

SAO-C5 


12 

3.986013 

Columbia 

Mercury 


13 

3.986013 

Columbia 

NAD -2 7 


Initial Datum Origin Shifts With Respect to NAD Origin (meters) (Reference 9) 


*Initial Datum 

A (meters) 

F 

AX 

AY 

AZ 

SAO-C5 

6378165.0 

298.25 

-31 

144 

181 

Mercury 

6378166.0 

298.30 

3 

111 

225 

NAD -2 7 

6378206.4 

294.9787 

0 

0 

0 


These initial coordinates are taken from References 9 and 10. The initial chords 
with respect to Columbia are also given for each datum. 

The Mercury datum ellipsoid and origin (not the individual site surveys) were 
derived from the NAD -27 datum ellipsoid and origin using all available data 
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including observations on Sputnik-1 and Vanguard (Reference 9). Thus, the 
Mercury datum coordinates contain both an origin translation and small sys- 
tematic individual station translations with respect to the NAD -27 datum co- 
ordinates. The size of the latter individual station translations is indicated by 
the slight differences in chord length from Columbia between the two datums. 

The SAO C-5 datum ellipsoid and origin were derived from SAO Baker Nunn 
camera observations on many satellites. In addition, the SAO camera site sur- 
veys were individually adjusted. The MOTS camera C-5 datum surveys were then 
obtained by a weighted interpolation of the SAO camera survey adjustments (Ref- 
erence 10). Therefore, the SAO C-5 datum coordinates also contain both an 
origin translation and small systematic individual station translations with respect 
to the NAD-27 datum coordinates. Again, the size of the individual station trans- 
lations is indicated by the slight differences in chords between the two datums. 

The chord differences between the SAO C-5 and the NAD-27 datums are generally 
several times larger than those between the Mercury and NAD-27 datums. 


4. RESULTS 

The results of the 13 solutions outlined in Table 1 are given in Tables 5 
through 17. For each solution the topocentric survey adjustments East, North, 
and Up (AE, An, AV) are given relative to the initial positions given in Tables 
2, 3, or 4. Also given in Tables 5 through 17 are the accuracy estimates of the 
adjustments based on an estimated accuracy of 1 arc second for the camera ob- 
servations and a priori estimates of 100 meters for the site coordinates for the 
12 adjusted stations. 

In each solution, the whole network, including the previously fixed origin 
station, is then shifted by the mean AE, An, AV adjustment for the network to 
provide a reasonable correction for the effects of an error in the origin station 
initial position. The resulting AE, AN, AV corrections are called the adjusted 
corrections. These adjusted corrections are also given in geocentric cartesian 
coordinates AX, AY, AZ. 

The chord changes, AC, are calculated as the difference between the adjusted 
chords and initial chords. The adjusted chord is the distance between the ad- 
justed station and adjusted origin positions. Also given are the proportional scale 
changes, S, due to these chord changes (S =^?). 

c 

Figures 5 through 17 illustrate the East (AE - AE) and North (AN - AN) 
vector components of the adjusted corrections given in Tables 5 - 17. The num- 
bers in parentheses on these figures are the station height adjustments 
(AV - AV). 
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The results of the 13 solutions are summarized below according to their 
effect on the parameters investigated, as indicated in Table 1. 

As a matter of interest, in some of the early SAO C-5 initial datum solutions, 
the initial height for Bermuda was incorrectly entered without a minus sign as 
+28.4 meters rather than -28.4 meters. These solutions recovered -55 to -68 
meters, of which 2 x 28.4 = 56.8 meters can be attributed to this input error, 
demonstrating the ability of the short arc techniques to successfully recover such 
errors. In Solutions 1 through 6, -56.8 meters were removed from the recovered 
Bermuda heights to avoid an erroneous effect on the chord calculations. 


(1) Effect of Error in^ 

The first four solutions (Tables and Figures 5-8) investigating the ef- 


fect of an error in p. are summarized below: 


Solution 

V- 

(xlO 1 4 m 3 /sec 2 ) 

Ap 

M2 

(X10 * ) 

AC 

(meters) 

S 

(xlO" 6 ) 

8 S 

(xio- 6 ) 

38 S 

1 

3.986032 

+4.77 

14.3 

8.3 

2.6 

7.8 

2 

3.986013 

0 

9.7 

5.7 

0 

0 

3 

3.985994 

-4.77 

7.0 

4.1 

-1.6 

-4.8 

4 

3.985956 

-14.30 

-2.8 

-1.3 

-7.0 

-21.0 


As is well known (Reference 11), a unit proportional change in p should pro- 
duce a 3-unit proportional change in scale. This comes from Kepler's Law: 


n 2 a 3 = M 


1 dM _ da _ 

so - — = oh 

3 p a 


(1) 


where it is assumed the change in scale of the orbits leads to a change in scale 
of the derived station positions. 


The first four solutions utilized different values of fi. The proportional 
changes inp relative to the JPL /z of solution 2 are given as The average 
chord adjustments and average proportional scale adjustments, obtained in each 

solution (Tables 5-8) are given under the headings AC and S where AC = ^EACj 

and S = Y2 SS i* The clian § es in s relative to the solution 2 adjustment are listed 

under 5S. These values are then multiplied by the factor 3 in order to compare 
directly with the values of4^. The results are plotted in Figure 3, and according 

*2 
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Figure 3. Effect of Error in /u 



ing to equation (1) should lie on a straight line of unit slope. The least squares 
fitted straight line of unit slope was calculated and is shown in Figure 3. This 
represents the theoretical behavior of equation (1), allowing for possible error 
in the reference solution (Solution 2). The largest residual with respect to this 
line represents an anomalous change in scale of only 1.9 x 10- 6 beyond the the- 
oretically expected change from equation (1). 

If the slope of the line is also allowed to adjust in the least squares fit, we 
obtain a slope of 1.5, as shown in Figure 3, and the largest remaining residual 
with respect to this new line represents a difference in scale of only 0.5 x 10 6 . 
This best fitting line of slope 1.5 is equivalent to a line of unit slope in a plot of 
28S rather than 35S versus ^ . The better fit of the steeper slope may reflect 
a departure from the region of linearity where equation (1) applies. In any case, 
these solutions verify that significant scale changes arise due to changes in jx. 
Thus, dynamic constraints do affect these short arc solutions, and it is impor- 
tant to choose the best available value for^u.. We have assumed the JPL value 
of 3.986013 x 10 14 m 3 /sec 2 to be the best available value. 


(2) Effect of Fixed Origin Choice 

Solutions 4, 5, and 6 (Tables and Figures 8 - 10) investigated the effect 
the choice of the fixed origin station has on the adjustment of the other stations. 
Origin stations were chosen so that one was near the network center, one east 
of this center, and one near the western limit of the network. 


The results of the three solutions are summarized below: 


Solution 

Origin 

AC 

(meters) 

S 

(x 10' 6 ) 

8 S 

(x 10~ 6 ) 

4 

Columbia 

-2.8 

-1.3 

0 

5 

Rosman 

-7.7 

-4.0 

-2.7 

6 

Mojave 

1.2 

0.9 

2.2 


As before, the average chord change is given as AC, the average scale 
change is S, the relative scale change with respect to solution 4 is 8S. As can 
be seen, choosing the origin towards the eastern or western edge of the network 
does result in small scale changes of a few parts per million relative to the 
solution with origin near the center of the network. The magnitude of these 
scale changes due to change in origin are about equivalent to the estimated ac- 
curacy (Reference 1) in the solutions for scale. From considerations of sym- 
metry and geometrical strength, it is probably better to utilize the central station 
as the fixed origin station. 
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(3) Effect of Error in Fixed Origin 

Solutions 7 through 11 (Tables and Figures 11 - 15) were made to de- 
termine the effect of an error in the position of the origin station on the other 
site adjustments. 

Solution 7 assumed the published survey for the Columbia station on the 
Mercury datum. Solution 8 was made with Columbia's X-component of its 
cartesian coordinates perturbed by -50 meters. This transforms into a primarily 
westward movement of the station site. Solutions 9 and 10 were made with the 
Y-component perturbed by -50 meters first and then by +50 meters. A decrease 
in the Y- component transforms into primarily an increase in station height. 

In solution 11, the Z-component was perturbed by -50 meters, which transforms 
into a predominantly southward shift and some decrease in station height. 

These solutions are summarized as follows: 




Tables 11 - 

15 

Tables 11A - 

15A 



. Summaries 

Summaries 

Solution 

Columbia 

Origin 

AC 

(meters) 

S 

(X10‘ 6 ) 

8S 

(xlO’ 6 ) 

AC 

(meters) 

S 

(x 10‘ 6 ) 

8S 

(x 10' 6 ) 

7 

x 0 , Y 0 , Zq 

6.3 

3.7 

0 

5.1 

2.8 

0 

8 

X 0 -50m, Y 0 , Z 0 

6.3 

3.6 

-0.1 

6.3 

3.4 

0.6 

9 

X 0 ,Y -50m, Z 0 

-19.8 

-8.3 

-12.0 

O 

L O 

1 

-3.1 

-5.9 


X o ,Y o +50m,Z o 

28.7 

16.2 

12.5 

15.8 

8.8 

6.0 


Xq ,Y 0 ,Z 0 -50m 


7.4 

3.7 

10.9 

6.2 

3.4 


The column labeled 8S is the change in scale for the various solutions rela- 
tive to solution 7. As can be seen, a 50-meter shift in the X-direction (West) 
of the origin station position causes only a -0.1 x 10" 6 change in scale relative 
to the unshifted position (solution 7). This type of shift can be absorbed by the 
orbital elements and by a lateral shift of the whole network and still maintain 
the observation directions in right ascension and declination. On the other hand, 
a 50-meter shift of Columbia in the ±Y-directions, corresponding primarily to 
^height shifts, leads to a ±12 x 10' 6 change in scale relative to solution 7. As 
noted, a decrease in Y of the origin station, corresponding to an increase in 
height, causes a corresponding decrease in scale or a shrinkage of the network 
and vice-versa. The 50-meter shift in the Z-direetion is predominantly south, 
but also significantly reduces the height component leading to a scale expansion 
of 3.7 x 10 6 . 

This situation is illustrated in Figure 4, where the observation directions 
to the satellite, S, are shown from the origin station, O, and from the observation 
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Figure 4. Change in Scale Due to Adjustment Method. 

stations A and B. Let it be assumed that the total set of observation directions 
from all stations in the net has the geometrical strength to maintain the config- 
uration of adjusted station positions relative to the origin station invariant under 
translations of the origin station. If this assumption is valid, a translation of the 
origin station should not change network scale. Then, it follows that an increase 
in height of the origin station from O to O, , moves the derived satellite position 
and station positions by the same translation from S to S 1 , A to A 1 and from 
BtoB, with no change in_scale. However, removing the mean biases in local 
topocentric shifts (AE, AN, AV) from the network station position corrections 
to obtain the adjusted corrections, as described earlier, tends to translate O, 
back to O, A, to A 2 , and B, to B 2 . This results, as shown in Figure 4, in a net 
decrease in scale, which is about the right magnitude to explain our results. 
Thus, if the above assumption is correct, the recovered scale changes are due to 
our shortcut technique of adjusting the derived corrections by AE, AN, AV. 

We should have converted each station local topocentric shift AE, AN, AV. 
to AX, AY, AZ in_a common cartesian system and used these to obtain 
mean biases AX, AY, AZ with which to adjust the corrections. 

The latter technique was applied to our solutions 7 - 11. The results are 
detailed in Tables 11A - 15A and summarized above. These solutions in general 
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contain about half the changes in scale due to AY origin station shifts as were 
previously obtained. However, some scale change still does occur and there- 
fore the assumption of scale invariance under origin station translations is not 
entirely valid. 


The conclusion drawn from these studies is that an error in origin station 
initial height relative to the initial surveys of the observation stations leads to 
scale changes on the order of -1.5 x 10 4 * 6 per +10 meter error in origin station 
height, assuming the station corrections are adjusted by Ax, AY, AZ, 
rather than by AE, AN, AV. 

height, 

= X = -0.15 x 10' 6 per meter). 

^ ' ■’ AH o - AY o cos <P 0 - 50 cos 39° 

The scale changes due to origin station height errors are approximately 
doubled by adjusting the corrections by AE, AN, Av rather than by AX, AY, 
A~Z. This indicates that for regular solutions, such as 2, 12, and 13, where the 
origin station average height adjustments AY, are only 3 to 4 meters, the un- 
wanted effect on network scale due to our method of adjustment is only about 
-0.5 x 10~ 6 . 


It is also of interest to note the ability of the solutions to recover the 50 
meter AX, AY, AZ errors injected into the origin station positions. The results 
are summarized below where the corrections are normalized relative to solution 
7 in order to remove the contributions of the data from those of the origin shifts. 
It can be seen that this recovery is good to 1.7 meters or better. 




Avg. Network Corrections 
(meters) 

Adjusted Avg. Network 
Corrections Relative to 
Solution 7 (meters) 

Solution 

Columbia 

Origin 

AX 

AY 

AZ 

AX 

AY 

AZ 

7 

X 0 , y q , Z 0 

4.2 

6.9 

2.2 

0 

0 

0 

8 

X 0 - 50m 

-44.6 

7.2 

2.8 

-48.8 

0.3 

0.6 

9 

Y 0 - 50m 

0.0 

-41.5 

7.6 

-4.1 

-48.4 

5.4 


Y 0 50m 

6.0 

55.4 

1.0 

1.9 

48.5 

-32 


Z 0 - 50m 

6.2 

6.3 

-48.3 

2.0 

-0.6 

-50.5 


(4) Effect of Choice of Initial Datum 

The last effect investigated was the choice of datum to which the initial 

station positions are referred. A major effect in using the different initial da- 
tums is the translation of all the initial station positions (not just the origin sta- 
tion) by the amounts AX, AY, AZ, given in Table 1 and repeated below. This is 
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equivalent to translating the Earth's Center of Mass by these values. The re- 
sults are given in solutions 2, 12, and 13, which are summarized below: 




Initial Datum Origin 
Shifts From NAD (meters) 

Results of Solutions 

Solution 

Initial Datum 

AX 

AY 

AZ 

■ 


ipplpl 

mm 

SAO-C5 

-31 

144 

181 

9.7 

5.7 

+2.4 


Mercury 

3 

111 

225 

2.8 

1.6 

-1.7 

13 

NAD -2 7 

0 

0 

0 

4.9 

3.3 

0 


The change in scale of 4.1 x 10' 6 obtained between solution 12 with the 
Mercury datum and solution 2 with the SAO-C5 datum is supported by a similar 
change of 4.6 x 10 6 observed between solutions 7 (Mercury datum) and 1 (SAO- 
C5 datum) while using a different value of/u.. 

The scale changes observed between solutions 2, 12, and 13, due primarily 
to the given center-of-mass translations, are about the same net magnitude as 
might be expected from the same translations of the fixed origin station, as 
described in the previous section. About 1/3 of the observed average change in 
scale is due to the small systematic differences in the Cartesian coordinates 
for the initial station positions. 

The SAO-C5 dynamically improved origin station initial position relative to 
the Earth Center of Mass is taken to be the best starting value for these short 
arc solutions. 


5. CONCLUSIONS 

1. The short arc solutions are affected in scale by the value chosen for 
p, more or less as anticipated by equation (1). The decrease in /u. from 
the Mercury value in solution 1 to the JPL value in solution 2, a pro- 
portional change of -4.77 x 10‘ 6 , causes a change in network scale of 
-2.6 x 10' 6 . Therefore, the best available value for/u should be used. 

2. The solutions are also affected by 2 to 3 x 10' 6 in scale by choosing 
the reference or origin station near the extremes of the network rather 
than near the more favorable central location. 

3. An error in height of the central origin station relative to the other 
station initial surveys causes a scale change of about -1.5 x 10 
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per + 10 meter height error. Therefore, the best available origin sta- 
tion local survey should be used. 

4. An error in origin station position relative to the Earth Center-of-Mass 
causes about the same scale change as the error relative to the local 
survey. Therefore, the best available Earth Center of Mass System 
should be used. 

5. In summary, solution 2 utilizes the best available value for/z,, a cen- 
trally located origin station, the best available local survey for the 
origin station, and a dynamically improved value for the origin station 
position relative to the Earth's Center of Mass, and is therefore the 
best solution in this study. It is recognized that more recent and prob- 
ably better dynamic solutions for the Earth's Center of Mass have be- 
come available since this study was initiated. 
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AVERAGE -1.9 -13.7 33.6 0.1 -41.5 7.6 -5.6 























Solution 10, Without Adjusted Corrections 
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AVERAGE 11.1 24.5 -40.5 6.1 55.4 -1.0 15.8 













Solution 11, Without Adjusted Corrections 
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Figure 6. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 2 
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- 20 Meters 

Figure 7. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 3 
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Figure 8. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 4 
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Figure 9. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 5 



Solution 6 



42 


Figure 10. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 6 
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Figure 11. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 7 
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Figure 12. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 8 
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Figure 13. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 9 
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Figure 14. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 10 
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Figure 15. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 11 



Solution 12 



48 


Figure 16. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 12 
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Figure 17. Horizontal Displacements Between Initial Survey and Satellite Survey — Solution 13 



